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Abstract 

We  begin  a  study  of  the  dynamics  of  ice  motion  during  river  breakup  by 
formulating  a  kinematic  model .  Ice  continuity  equations  ore  applied  to  relate  the 
speeds  of  a  breaking  front,  convergence  front,  stoppage  front,  and  release  front 
with  the  ice  discharge  and  volume  per  unit  surface  area  (unit  volume)  on  either 
side  of  each  front.  Ice  velocity  data  are  obtained  from  measurements  with  time 
made  during  a  dynamic  breakup  at  a  pair  of  sites  bounding  a  reach  of  the 
Connecticut  River.  We  simulate  the  ice  and  front  motion  through  time  for  this 
reach  using  the  kinematic  model  with  the  assumptions  that  accumulation 
thickness  and  porosity  are  uniform,  and  that  changes  in  the  ice  conditions  and 
motion  occur  only  at  a  front.  Contrary  to  the  basic  assumption  of  static  jam 
formation,  we  find  that  the  accumulation  develops  while  the  ice  is  moving,  and 
that  jam  formation  merely  represents  the  arrest  of  the  motion. 


Cover:  (Top)  Setup  of  targets  on  a  known  grid  for  video  ice  velocity  measure¬ 
ment  at  site  Don  the  Connecticut  River  at  Cornish,  New  Hampshire.  Ice 
breakup  occurred  one  week  later.  (Bottom)  Looking  across  the  river 
during  brash  ice  motion  at  about  1  m/s.  (Photos  by  M.  Ferrick.) 


For  conversion  of  SI  metric  units  to  U.S./British  customary  units  of  measurement 
consult  ASTM  Standard  E380-89a,  Standard  Practice  for  Use  of  the  International 
System  of  Units,  published  by  the  American  Society  for  Testing  and  Materials, 
1916  Race  St.,  Philadelphia,  Pa.  19T03. 
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Kinematic  Model  of  River  Ice  Motion 
During  Dynamic  Breakup 

MICHAEL  G.  FERRICK,  PATRICIA  B.  WEYRICK  AND  DAVID  F.  NELSON 


INTRODUCTION 

The  process  of  dynamic  ice  breakup  in  a  river  is 
complex  because  of  the  interdependence  of  the  wa¬ 
ter  flow  and  the  ice  motion,  and  the  rapid  changes 
with  time  in  both  of  these  processes.  The  unsteady 
flow  associated  with  dynamic  breakup  has  been  con¬ 
sidered  by  several  authors  (e.g.,  Beltaos  and 
Krishnappan  1982,  Billfalk  1982,  Doyle  and  Andres 
1979,  Ferrick  and  Mulherin  1989,  Henderson  and 
Gerard  1981,  Prowse  et  al.  1986,  Williamson  1989). 
However,  in  the  static  analyses  of  breakup  ice  jan. 
formation  of  Beltaos  (1988),  Pariset  et  al.  (1966)  and 
Uzuner  and  Kennedy  (1976),  die  flow  and  ice  dynam¬ 
ics  are  neglected  and  jam  formation  is  assumed  to 
be  independent  of  other  breakup  processes.  The  static 
force  balance  requires  the  ice  to  thicken  in  place  by 
shoving  and  collapse  until  an  equilibrium  is  reached 
between  the  internal  strength  of  the  accumulation 
and  the  steady  hydraulic  forces.  Gerard  and  Flato 
(1988)  emphasized  the  need  for  understanding  the 
hydraulics  and  mechanics  of  moving  ice  accumula¬ 
tions  and  jam  formation  in  the  context  of  dynamic 
breakup.  Recently  Shen  and  Chen  (1992)  developed 
a  dynamic  model  of  flow  and  ice  motion  applied  to 
the  ffeezeup  processes  of  a  river,  and  Guo  (1991)  pro¬ 
posed  a  dynamic  model  of  ice  breakup. 

A  first  step  in  describing  many  river  ice  and  flow 
processes  is  the  conservation  of  mass.  Lighthill  and 
Whitham  (1955)  used  continuity  to  develop  a  kine¬ 
matic  wave  model  of  traffic  flow  on  crowded  roads. 
When  congested  traffic  ahead  is  slower  than  that  be¬ 
hind,  these  waves  coalesce  into  a  "shock  wave"  or 
front  that  causes  rapid  speed  reduction  of  vehicles 
upon  arrival.  Risebro  and  Tveito  (1992)  presented  a 
front  tracking  method  for  solving  hyperbolic  conser¬ 
vation  laws.  The  solution  is  represented  by  an  arbi¬ 
trary  number  of  piecewise  constant  states  separated 
by  discontinuities  or  fronts.  These  fronts  are  shock 


or  rarefaction  waves,  whose  motion  is  tracked  in  the 
x-t  plane  until  they  interact. 

In  this  report  we  obtain  simultaneous  ice  velocity 
records  with  time  at  a  pair  of  sites  bounding  a  reach 
of  the  Connecticut  River  from  fixed  video  cameras 
and  orthogonal  target  grids  (Ferrick  et  al.  1992).  The 
ice  motion  at  the  sites  differed,  raising  questions  that 
include  how  these  views  of  the  motion  are  related, 
and  whether  the  data  obtained  at  two  points  can  be 
used  to  understand  the  ice  dynamics  of  the  reach. 
An  analogy  is  proposed  between  the  ice  motion  dur¬ 
ing  dynamic  breakup  and  the  flow  of  traffic  on 
crowded  roads.  We  then  develop  a  kinematic  model 
of  river  ice  motion  from  the  continuity  equation  and 
identify  several  fronts  in  the  ice.  The  fronts  are 
tracked  in  the  x-t  plane  to  obtain  a  hypothetical  event 
sequence  and  to  describe  the  ice  motion  in  the  reach 
from  the  data.  The  resulting  ice  accumulation  forms 
while  the  breakup  is  in  progress.  The  dynamic  track¬ 
ing  of  ice  particles  and  fronts  with  jamming  as  an 
integral  part  of  breakup  replaces  the  usual  assump¬ 
tion  of  static  jam  formation. 

KINEMATIC  MODEL 
OF  ICE  MOTION 

To  construct  our  model,  we  will  visualize  the  ex¬ 
istence  of  different  "fronts"  that  move  through  the 
ice  field  in  response  to  changing  conditions.  The 
breaking  front  travels  downstream  and  separates  the 
stationary  ice  ahead  from  the  moving  ice  behind.  The 
convergence  front  moves  with  the  upstream  limit  of 
an  ice  accumulation,  delineating  it  from  thinner  ice 
upstream.  Analogies  for  these  fronts  occur  at  both 
ends  of  a  congested  traffic  accumulation  moving  at 
a  reduced  speed.  Additional  vehicles  regularly  added 
to  the  front  of  the  accumulation  cause  the  breaking 
front  to  travel  faster  than  die  group.  Behind  the  group 


at  the  convergence  front,  faster  moving  vehicles  slow 
to  the  group  speed  and  become  part  of  the  accumu¬ 
lation.  The  stoppage front  travels  upstream  following 
the  arrest  of  ice  motion,  and  separates  the  moving 
ice  upstream  from  ice  at  rest  downstream.  The  re¬ 
lease front  moves  upstream  through  the  ice  field,  ini¬ 
tiating  motion  in  response  to  an  ice  release  down¬ 
stream.  A  red  light  interrupts  a  traffic  flow,  causing 
vehicles  to  stop  in  response  to  those  ahead.  The  stop¬ 
page  front  progresses  "upstream"  at  a  rate  depen¬ 
dent  on  vehicle  separation.  Following  a  green  light, 
vehicle  motion  is  initiated  in  response  to  that  ahead, 
and  this  release  front  also  travels  "upstream." 

Ice  continuity  equation 

The  conservation  of  mass  can  be  written  for  the 
ice  near  each  front  that  relates  the  speed  of  the  front 
to  the  change  in  ice  thickness  and  accumulation  po¬ 
rosity  across  the  front.  The  ice  discharge  on  either 
side  of  each  front  can  also  be  obtained.  Ice  continu¬ 
ity  equates  the  time  rate  of  change  of  mass  in  an  ar¬ 
bitrary  control  volume  (CV)  to  the  net  efflux  of  mass 
across  the  control  surface.  We  will  use  the  CV  de¬ 


picted  in  Figure  1,  containing  and  moving  with  a 
front  at  speed  C  in  the  downstream  direction.  Ice  does 
not  accumulate  in  this  moving  CV  and  there  is  no 
time  rate  of  change  term  to  consider.  Ice  moving  at 
velocity  Vq  enters  the  control  volume  through  down¬ 
stream  control  surface  1,  and  an  equal  quantity  exits 
the  control  volume  through  upstream  surface  2  at 
Vq.  The  relative  velocities  between  the  ice  and  the 
upstream  and  downstream  control  surfaces  are  C  - 
Vj2  and  C  -  Va,  respectively.  The  areas  of  the  control 
surfaces  are  8t„  where  B  is  the  river  width,  and  f,  is 
the  thickness  of  the  ice  sheet  or  accumulation.  The 
ice  continuity  equation  then  states  the  flux  balance 
as 

Bi  (C- V,i)  hi  (l-eci)  = 

B2(C-Vi2)ti2{  l-e<2)  (1) 

where  ec  is  the  porosity  of  the  ice  accumulation.  We 
refer  to  the  ice  volume  per  unit  surface  area,  f,  (l-ec), 
as  the  unit  ice  volume.  The  general  form  of  eq  1  is 

(C-Vn)  =  R(C-Vi2)  (2) 


Moving 

Control  Volume 


Figure  1.  Control  volume  with  control  surfaces  1  and  2  that  either  moves  downstream  or 
upstream  with  the  corresponding  front.  An  abrupt  parameter  transition  at  a  front  between 
constant  states  on  each  side  of  the  CV  is  shown  for  ice  thickness. 
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Bl  t|i  (1-flci)  ^ is 

Figure  2.  Dimensionless  breaking  front  speed  as  a  function  of  (left)  relntive  river  width 
times  unit  ice  volume  on  either  side  of  the  front,  and  (right)  relative  width  times  ice  thickness 
for  incoming  rheet  and  selected  values  of  accumulation  porosity. 

where  R  represents  the  dimensionless  ratio  of  effec-  =  ^b.  (4) 

tive  control  surface  areas  as  V\2  Rb  ~  1 

B  t  (l-e  )  with  R  at  the  breaking  front,  Rb  >  1.  If  the  breaking 

R  =  -  - - R- .  (3)  front  encounters  an  unbroken  ice  sheet  of  thickness 

fii  ( 1  -ec\)  instead  of  ice  rubble,  =  fis  and  if  there  is  negli¬ 

gible  open  water  area  in  the  sheet,  the  incoming  sheet 
In  the  development  of  eq  2  the  CV  was  assumed  to  "porosity"  ecl  =  0. 

move  in  the  downstream  direction.  This  sign  con-  The  movement  of  a  breaking  front  causes  ice  con¬ 
vention  yields  negative  speed  -C  for  fronts  moving  vergence,  which  results  in  a  decreased  surface  area 

upstream.  occupied  by  the  ice.  The  relationship  between 

dimensionless  breaking  front  speed  and  relative  river 
Breaking  front  width  times  unit  ice  volume  given  by  eq  4  is  depicted 

We  will  first  consider  a  reach  of  broken  ice  that  in  Figure  2.  Significant  changes  in  river  width  do  not 

has  jammed  upstream  of  an  intact  ice  cover.  A  flow  usually  occur  in  a  longitudinal  distance  of  less  than 

surge  traveling  downstream  arrives  at  the  accumu-  one  river  width.  When  the  ice  passing  through  the 
lation,  and  at  some  point  ice  motion  is  initiated.  control  volume  thickens  substantially  or  the  poros- 
Downstream  ice  movement  continually  adds  mate-  ity  of  the  accumulation  decreases  greatly,  the  break- 
rial  to  the  front  of  the  moving  pack,  and  the  speed  of  ing  front  speed  is  only  marginally  higher  than  the 

the  breaking  front  Cb  is  always  greater  than  the  local  ice  velocity.  However,  as  the  ratio  of  unit  ice  volumes 

ice  velocity  Va.  With  appropriate  substitutions,  the  across  the  front  approaches  1,  the  breaking  front  ce- 
dimensionless  breaking  front  speed  can  be  obtained  lerity  becomes  large.  Figure  2  also  presents  the  di- 

from  eq  2  as  mensionless  breaking  front  speed  of  advance  into  an 
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Figure  3.  Dimensionless  ice  conver¬ 
gence  front  speed  as  a  function  of  rela¬ 
tive  river  width  times  unit  ice  volume 
on  either  side  of  the  front  fora  range  of 
V,  j/Vl2.  Negative  or  upstream  speeds  are 
distinguished  from  positive  downstream 
speeds. 

intact,  stationary  ice  sheet,  as  a  function  of  the  accu¬ 
mulation  to  sheet  thickness  ratio  for  a  range  of 
accumulation  porosities.  The  thickness  ratio  at  which 
the  breaking  front  speed  becomes  very  large  depends 
on  the  porosity  of  the  accumulation.  For  larger  thick¬ 
ness  ratios  this  speed  is  small  for  all  porosities. 

Convergence  front 

In  the  ice  accumulation  immediately  upstream  of 
the  breaking  front  the  total  resistance  to  motion  is 
relatively  high.  Farther  upstream,  the  motion  resis¬ 
tance  is  lower  and  higher  ice  velocities  are  typical. 
This  velocity  difference  causes  ice  convergence  at  the 
upstream  limit  of  a  moving  accumulation.  For  sim¬ 


plicity  we  visualize  this  accumulation  as  a  rigid  body 
with  a  uniform  velocity.  All  ice  convergence  occurs 
at  die  convergence  front,  located  at  the  upstream  limit 
of  the  accumulation.  Consider  a  CV  that  moves  with 
the  convergence  front  at  the  head  of  the  accumula¬ 
tion  at  speed  Cc.  The  control  surface  designations  are 
those  of  Figure  1,  and  the  dimensionless  speed  of  the 
ice  convergence  front  is  obtained  from  eq  2  as 

^i-Rc 

=  Yi2 -  (5) 

Vi2  1-RC 

where  is  the  accumulation  velocity,  is  the  ve¬ 
locity  of  the  ice  immediately  upstream,  Rc  <  1  indi¬ 
cates  convergence  of  thinner  or  more  porous  ice  from 
upstream  of  the  accumulation,  and  Cc  is  positive  for 
downstream  front  motion  or  negative  for  fronts  mov¬ 
ing  upstream.  The  dimensionless  convergence  front 
speed  is  given  in  Figure  3  as  a  function  of  Rc  for  a 
range  of  ice  velocity  ratios.  Ice  accumulations  up¬ 
stream  of  a  relatively  slow  moving  breaking  front 
would  generally  have  velocity  ratios  of  less  than  1. 
Then,  as  the  speed  of  the  accumulation  decreases,  the 
convergence  front  speed  also  decreases  or  becomes 
more  negative  for  all  Rc  values.  This  front  speed  also 
decreases  or  becomes  more  negative  as  Rc  increases 
for  all  velocity  ratios.  However,  divergence  and 
shortening  at  the  upstream  limit  of  the  accumulation 
are  indicated  by  dimensionless  front  speeds  greater 
than  1,  and  the  front  speed  increases  with  both  Rc 
and  the  velocity  ratio. 

The  growth  of  an  ice  accumulation  during  a  pe¬ 
riod  of  steady  motion  is  shown  on  the  x-t  plane  in 
Figure  4.  Paths  of  front  motion  and  selected  ice  mo¬ 
tion  are  displayed  in  time  and  space  on  this  plane. 
Stationary  ice  initially  downstream  of  the  breaking 
front  is  incorporated  into  the  accumulation  as  the 
front  passes.  Ice  initially  upstream  of  the  convergence 
front  overtakes  it  and  is  also  incorporated  into  the 
accumulation.  The  total  length  of  the  accumulation 
is  the  x-distance  between  these  fronts.  Both  the  break¬ 
ing  and  convergence  fronts  deflect  the  path  lines  trac¬ 
ing  the  ice  motion.  The  parallel  path  lines  between 
the  fronts  reflect  assumed  rigid  body  motion  of  the 
accumulation. 

Stoppage  front 

A  continuity  equation  can  be  written  for  the  stop¬ 
page  of  ice  motion  that  progresses  upstream  follow¬ 
ing  the  arrest  of  a  breaking  front  or  of  the  moving 
ice  somewhere  behind  this  front.  In  this  case  the  CV 
depicted  in  Figure  1  moves  upstream  at  speed  Cs 
with  the  stoppage  of  ice  motion.  Moving  ice  enters 
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Figure  4.  Steady  motion  and 
growth  of  an  ice  accumulation 
(shaded)  depicted  on  the  x-t 
plane.  The  spatial  and  temporal 
separation  of  ice  particles  is  re¬ 
duced  upon  entering  the  accumu¬ 
lation  and  constant  thereafter. 


Figure  5.  Dimensionless  ice  stoppage  front  speed  and  di¬ 
mensionless  ice  release front  speed  (upstream)  as  a  function 
of  relative  river  width  times  unit  ice  volume  on  either  side  of 
the  front. 


through  the  upstream  control  surface  and  exits 
through  the  downstream  surface  as  stationary  ice. 
Solving  eq  2  for  the  dimensionless  speed  of  the  ice 
stoppage  front  yields 


^  2  1-Ks 


with  Rg  <  1,  indicating  the  unit  ice  volume  down¬ 
stream  of  the  front  is  greater  than  that  upstream. 
Therefore,  packing  during  an  ice  stoppage  is  limited 
by  the  downstream  unit  ice  volume  that  cannot  be 
exceeded.  Low  stoppage  front  speeds  in  Figure  5  re¬ 
flect  significant  packing  during  ice  arrest,  but  as  Rj 
approaches  1,  the  front  progresses  upstream  very 
rapidly.  The  breaking  and  stoppage  fronts  are  op¬ 
posites,  and  their  continuity  equations  have  the  same 
form.  The  distinction  is  the  inverse  relationship  be¬ 
tween  the  unit  ice  volume  across  each  front. 

A  stoppage  front  is  the  limiting  case  of  a  conver¬ 
gence  front  with  a  velocity  ratio  of  zero.  As  the  speed 
of  an  ice  accumulation  approaches  zero,  Cc  becomes 
negative,  the  convergence  front  becomes  a  steppage 
front,  and  eq  5  reduces  to  eq  6.  An  ice  stoppage  is 
illustrated  on  the  x-t  plane  in  Figure  6.  The  initial 
stoppage  occurs  upstream  of  the  breaking  front,  and 
the  stoppage  front  travels  upstream  instantaneous¬ 
ly  through  the  rigid  body  accumulation.  The  stop¬ 
page  and  convergence  fronts  merge  into  a  single  front 
at  the  upstream  limit  of  the  accumulation.  Then,  the 
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Figure  6.  Ice  stoppage  occurring  upstream  of  the  breaking  front  depicted  on  the 
x-t  plane.  The  path  lines  show  reduced  spatial  separation  of  ice  particles  upstream  of 
the  accumulation  following  passage  of  the  stoppage  front. 


stoppage  front  continues  upstream  at  speed  Cv 
bringing  that  ice  to  a  halt. 

Release  front 

We  can  also  write  a  continuity  equation  for  a  fail¬ 
ure  of  the  resistance  at  a  point  and  the  initiation  of 
ice  motion.  In  this  case  the  CV  shown  in  Figure  1 
moves  upstream  at  a  speed  C;  with  the  release  of  the 
ice  and  initial  motion.  Stationary  ice  enters  through 
the  upstream  control  surface,  and  exits  through  the 
downstream  surface  as  moving  ice.  The  dimension¬ 
less  speed  of  the  ice  release  front  is  obtained  from  eq 
2  as 


with  R;  >  1.  Ice  release  causes  divergence,  as  the  wa¬ 
ter  surface  area  available  to  an  ice  accumulation  is 
increased.  The  dimensionless  speed  of  the  ice  release 
front  is  depicted  in  Figure  5.  With  R,  =  1  ice  diver¬ 
gence  is  small,  approximating  rigid  body  motion,  and 
the  front  speed  is  very  large.  However,  as  R,  in¬ 
creases,  indicating  substantial  divergence,  the  dimen¬ 
sionless  release  front  speea  approaches  zero.  An  ice 
release  is  depicted  on  the  x-t  plane  in  Figure  7.  A 
breaking  front  travels  downstream  and  a  release 
front  travels  upstream  from  the  point  of  the  release. 


Both  fronts  initiate  ice  motion,  with  convergence  oc¬ 
curring  upstream  of  a  breaking  front  and  divergence 
occurring  downstream  of  a  release  front. 

Ice  discharge 

The  ice  discharge  at  point  p  Qp  can  be  written  as 

Qp  =  Bp  v\p  fip(l-^cp)  (8) 

Each  of  the  continuity  equations  associated  with  front 
motion  can  be  solved  to  obtain  the  ice  discharge  near 
the  front.  The  ice  discharge  upstream  of  th^  break¬ 
ing  front  Qf,  is  obtained  from  eq  4  and  8  as 

Qb  =  cb  Db  (9) 

where 

D  =  #2*12(1  -ea)-  Bifii(l  -ec\)  (10) 

and  Db  is  the  value  of  D  associated  with  the  break¬ 
ing  front.  The  ice  discharge  upstream  of  the  conver¬ 
gence  front  Qc  can  be  obtained  from  eq  5  and  8  as 

Qc~Qb  =  CcDc  (11) 

The  ice  discharge  upstream  of  the  stoppage  front  Qs 
can  be  obtained  from  eq  6  and  8  as 
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Figure  7.  Ice  release  at  a  point  initiating  breaking 
and  release  front  motion  depicted  on  the  x-t  plane. 
The  path  lines  show  increased  spatial  separation  of  ice 
particles  fallowing  passage  of  the  release  front. 


Qs  =  -QDS  (12) 

Finally,  the  ice  discharge  downstream  of  a  release 
front  Q,  is  obtained  from  eq  7  and  8  as 

Qi=CiDi  (13) 

1992  CONNECTICUT  RIVER 
BREAKUP  DATA 

Breakup  of  the  Connecticut  River  ice  cover  near 
Windsor,  Vermont,  occurred  on  11  March  1992.  Prior 
to  breakup  we  installed  an  orthogonal  grid  of  video 
targets  with  known  dimensions  at  two  sites,  sepa¬ 
rated  by  a  1600-m  reach.  The  upstream  site  was  des¬ 
ignated  site  U  and  the  downstream  site  was  D.  Tri¬ 
pod-mounted  video  cameras  were  positioned  at  the 
top  of  the  river  bank  to  record  the  motion  of  the  ice 
through  the  grid  at  each  location.  The  cameras  re¬ 
mained  focused  on  the  original  grid  location  for  the 
entire  period  of  motion,  and  audio  marks  on  each 
videotape  were  used  to  synchronize  the  times.  Each 
videotape  was  window  dubbed  with  continuous  on¬ 
screen  digital  time  accurate  to  ±0.03  s.  Points  on  the 
ice  to  be  used  for  velocity  measurement  were  cho¬ 
sen  by  their  contrast  from  the  surrounding  ice.  The 
grid  scale  divided  by  the  travel  time  provides  an  ice 
velocity  estimate  at  the  mean  time  of  the  interval.  The 
breakup  occurred  in  two  ice  movements,  separated 


by  about  3  hours  during  a  period  of  increasing  dis¬ 
charge.  In  this  report  we  will  analyze  the  initial  pe¬ 
riod  of  ice  motion  and  use  data  from  the  second  mo¬ 
tion  for  hypothesis  testing. 

Orthogonal  polynomial 
fits  to  the  data 

The  methods  described  by  Ferrick  et  al.  (1992) 
were  used  to  obtain  least-squares  orthogonal  poly¬ 
nomial  fits  to  the  ice  velocity  data.  The  polynomials 
are  used  to  quantify  the  variability  in  the  data  for 
comparison  with  expected  measurement  error  and 
to  verify  that  velocity  variations  across  the  channel 
can  be  neglected.  The  length  of  ice  passing  between 
specified  times  and  the  ice  acceleration  at  a  given  time 
are  obtained  by  integration  and  differentiation,  re¬ 
spectively,  of  the  corresponding  polynomial.  Accu¬ 
rate  estimates  of  ice  acceleration  cannot  be  obtained 
from  the  raw  data.  Polynomials  of  increasing  degree 
incorporate  additional  information  from  the  data  un¬ 
til  the  magnitude  of  the  maximum  error  ceases  to 
decrease  and  the  structure  of  the  error  is  random.  A 
polynomial  of  appropriate  order  can  be  found  by 
observing  both  the  absolute  magnitude  and  ran¬ 
domness  of  the  error.  The  7th-  and  12th-order  poly¬ 
nomials  are  presented  together  with  the  data  from 
site  D  in  Figure  8.  The  7th-order  polynomial  repre¬ 
sents  the  highest  order  for  which  a  composite  poly¬ 
nomial  equation  is  available.  The  errors  between 
these  polynomials  and  the  data  are  also  given  in  Fig- 
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Ice  Velocity  (m/s)  Error  (m/s) 


Figure  8.  Ice  velocity  data,  best  least- 
squares  7th  and  12th  order  polynomial 
fits,  and  velocity  error  between  the 
polynomials  and  the  data  for  site  D.  The 
predominant  floe  size  during  this  ice  mo¬ 
tion  was  on  the  order  of  100m. 


Table  1:  Errors  associated  with  polynomial  fits  to  the  ice  velocity  data. 


Site 

Polynomial 

degree 

Maximum 
polynomial 
ice  velocity 
(m/s) 

Maximum 

measurement 

error 

(m/s) 

Maximum 
error  (poly¬ 
nomial-data) 
(m/s) 

Mean  poly¬ 
nomial  ice 
velocity 
(m/s) 

Measurement 

error  at  mean 
velocity 
(m/s) 

RMS  error 
(polynomial- 
data) 
(m/s) 

Upstream 

7 

1.36 

0.109 

0.113 

0.89 

0.039 

0-420s 

0.065 

11 

134 

0.107 

0.095 

0.89 

0.034 

Upstream 

7 

2.37 

0322 

0.429 

134 

0.083 

420-970s 

0.097 

15 

2.34 

0318 

0.194 

134 

0.058 

Downstream 

7 

0.93 

0.069 

0.165 

0.66 

0.064 

0-420s 

0.047 

12 

1.02 

0.077 

0.089 

0.66 

0.028 
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Figure  9.  Ice  velocity  data,  best  least- 
squares  7th  and  11th  order  polynomial 
fits,  and  velocity  error  between  the  polyno¬ 
mials  and  the  data  for  the  first  part  of  the 
motion  at  site  U. 


ure  8.  The  7th-order  error  has  structure,  and  is  larger 
in  magnitude  than  the  12th-order  error.  Figure  9  pre¬ 
sents  the  data,  polynomials,  and  errors  for  the  first 
part  of  the  longer  record  at  site  U,  and  Figure  10  gives 
comparable  information  for  the  second  part  of  this 
record.  The  order  of  the  polynomial  providing  the 
best  least-squares  fit  depends  on  the  length  and  com¬ 
plexity  of  the  data  record.  The  best  polynomials  for 
each  time  segment  were  obtained  by  assuming  a 
match  point  at  a  common  time.  We  present  the  com¬ 
plete  data  record  for  site  U  and  a  continuous  poly¬ 
nomial  representation  of  ice  velocity  in  Figure  11. 

Maximum  measurement  errors  in  obtaining  data 
from  the  videotape  were  1 .0  m  in  length  and  0.2  s  in 
elapsed  time.  These  estimates  were  obtained  from 
known  target  size  and  width  of  grid  lines  superposed 


on  the  videotape,  and  the  precision  of  viewing  a  tar¬ 
get-line  intersection.  When  these  errors  are  combined 
unfavorably  the  error  in  the  computed  ice  velocity 
increases  with  the  ice  velocity.  A  maximum  error  of 
0.08  m/s  at  an  ice  velocity  of  1  m/s  increases  to  0.23 
m/s  at  a  velocity  of  2.5  m/ s.  Ice  velocity  errors  be¬ 
tween  the  various  polynomials  and  the  data  are  com¬ 
pared  with  measurement  errors  in  Table  1 .  The  maxi¬ 
mum  polynomial-data  errors  are  slightly  less  than 
the  maximum  measurement  errors  for  the  upstream 
site,  and  slightly  greater  for  the  downstream  site.  The 
RMS  polynomial-data  errors  average  only  57%  of  the 
maximum  measurement  error  at  the  mean  velocity 
and  32%  of  the  overall  maximum  measurement  er¬ 
ror.  Therefore,  most  of  the  error  in  the  ice  velocity 
polynomials  (Fig.  8, 9, 10)  can  be  attributed  to  mea- 
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Ice  Velocity  (m/s)  Ice  Velocity  (m/s)  Error  (m/s) 


Figure  10.  Ice  velocity  data,  best  least- 
squares  7th  and  15th  order  polynomial 
fits,  and  velocity  error  betiveen  the  poly¬ 
nomials  and  the  data  for  the  second  part  of 
the  motion  at  site  U.  A  piecewise  linear  ice 
velocity  for  the  motion  calculations  is  also 
shown. 


Figure  11.  Complete  ice  velocity  data 
record  for  site  U  and  best  11th  and  15th 
order  polynomials  with  a  match  point 
connecting  them.  The  predominant  floe 
size  is  indicted  as  a  function  of  time. 
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surement  error,  and  these  polynomials  represent  the 
ice  velocity  as  well  as  the  data.  There  is  no  evidence 
of  significant  velocity  differences  across  the  channel. 

Initial  breaking  front  motion 

The  river  discharge  prior  to  ice  motion  was  980 
m3/s  and  gradually  increased.  The  initial  ice  move¬ 
ment  was  the  result  of  a  support  breakup  of  the  river 
reach.  The  bank  support  failed  and  die  ice  sheet 
moved  as  a  large  plate  with  minimal  convergence 
and  rubble  formation.  Other  characteristics  of  sup¬ 
port  breakup  are  given  by  Ferrick  and  Mulherin 
(1989).  The  ice  near  the  banks  at  the  downstream  site 
foiled  almost  continuously  during  the  ice  motion.  The 
ice  acceleration  calculated  for  each  site  using  the  ve¬ 
locity  polynomial  was  within  3%  of  0.0266  m/s2  fol¬ 
lowing  breaking  front  arrival.  At  50  s  the  ice  veloci¬ 
ties  were  equal  at  the  sites,  an  indicator  of  nearly  rigid 
plate  initial  motion.  The  time  lag  between  initial  ice 
motion  at  the  two  camera  sites  was  only  2  s,  indicat¬ 
ing  a  breaking  front  speed  of  800  m/s  in  the  down¬ 
stream  direction.  With  QyV,  very  large,  the  moving 
control  volume  analysis  indicates  that  the  ice  sheet 
was  unchanged  as  the  breaking  front  passed  through 
the  reach  except  that  motion  was  initiated.  This  re¬ 
sult  concurs  with  the  analysis  of  Ferrick  et  al.  (1992) 
that  ice  upstream  of  a  rapidly  moving  breaking  front 
moves  almost  as  a  rigid  body.  The  fast  breaking  front 
speed  also  supports  the  modeling  concepts  of  force 
balance  per  unit  length  and  instantaneous  support 
breakup  of  a  reach  (Ferrick  and  Mulherin  1989). 

Water  fountaining  above  the  level  of  the  ice,  pos¬ 
sibly  indicating  a  pressure  flow,  was  not  observed 
near  either  site.  The  ice  cover  was  cracked  continu¬ 
ously  along  both  banks,  and  moved  readily  in  the 
vertical  direction  in  response  to  changing  water  lev¬ 
els.  Therefore,  the  rapid  front  motion  was  not  caused 
by  pressure  waves  in  the  flow.  Bourbie  et  al.  (1987) 
present  pressure  wave  speed  Vp  for  an  isotropic 
material  as 

(i4) 

Lp(l+v)  (1-2v)J 

where  £  is  the  elastic  modulus,  v  is  Poisson's  ratio, 
and  p  is  density.  The  pressure  wave  speed  in  fresh¬ 
water  ice  at  its  melting  point  varies  between  3700 
m/s  for  intact  ice  and  2000  m/s  for  deteriorated  ice. 
The  most  likely  cause  of  the  rapid  front  motion  was 
pressure  transmitted  downstream  through  the  ice, 
interrupted  and  slowed  occasionally  by  ice  failure. 


Overall  ice  mass  balance 

The  length  of  ice  sheet  that  passed  each  site  dur¬ 
ing  this  ice  motion  can  be  found  by  integrating  the 
ice  velocity  polynomials.  The  quantity  of  ice  that 
accumulated  in  the  reach  Aq  between  times  tQ  and  t} 
can  then  be  obtained  as 


Aq=Buf  Vmtiu{l  -ecJdt  -  Bd 
'to 


- 

Vidtid(l-ecd)dt  (15) 

■to 

where  subscripts  u  and  d  indicate  sites  U  and  D,  re¬ 
spectively.  The  river  width  at  site  U,  Bu  is  140  m,  and 
the  width  at  site  D,  Bd  is  135  m.  The  length  of  ice 
rubble  in  the  reach  Lr  is  measured  downstream  from 
site  U  to  the  breaking  front.  The  final  rubble  accu¬ 
mulation  contains  the  net  ice  accumulation  in  the 
reach  during  the  motion  and  the  additional  rubble 
that  was  generated  in  the  reach  at  the  breaking  front 
and  the  convergence  front 

Mir(l-eCr)  =  Mis  +  /lq/B  fl®) 

where  the  average  width  of  the  river  in  the  accumu¬ 
lation  reach,  B  is  190  m  and  the  subscript  r  indicates 
parameters  of  the  rubble.  Solving  eq  16  for  the  length 
of  the  rubble  reach  we  obtain 


Li-  V* 

fir(l  —  ecr)  ~  t\s 


(17) 


Then,  the  volume  of  ice  rubble  in  the  final  accumu¬ 
lation  Av  is  obtained  as 

Av=Btir(l-ecr)Lr  (18) 

Beltaos  (1988)  and  Calkins  (1978)  reported  exten¬ 
sive  ice  jam  thickness  measurements.  The  thickness 
was  generally  greatest  near  the  toe  and  decreased 
with  distance  upstream.  However,  in  all  cases  the 
mean  cross-sectional  thicknesses  were  nearly  con¬ 
stant  over  much  of  the  jam  length.  The  mean  thick¬ 
ness  of  these  jams  ranged  from  2  to  5  times  the  origi¬ 
nal  sheet  thickness.  The  ice  in  the  reach  between  the 
camera  sites  is  idealized  as  having  two  components: 
ice  plates  with  a  unit  volume  comparable  to  the  origi- 
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Table  2:  Parameters  for  a  range  of  accumulation  unit  volume. 


(m) 

Wo 

(m) 

4 

(m) 

Percent  of  reach 
as  rubble 

A„ 

(m3  x  ICT3) 

Rb 

<W* 

0.75 

120 

1100 

69 

157 

15 

3.0 

0.83 

90 

840 

53 

133 

1.67 

25 

1.00 

60 

570 

36 

108 

2.0 

2.0 

1.10 

50 

480 

30 

100 

2  2 

1.83 

1.20 

43 

410 

26 

94 

2.4 

1.71 

1.25 

40 

390 

24 

93 

25 

1.67 

1.30 

38 

370 

23 

91 

2.6 

1.63 

1.40 

33 

330 

21 

88 

2.8 

1.56 

1.50 

30 

300 

19 

86 

3.0 

15 

nal  sheet,  and  rubble  with  a  greater  unit  volume.  For 
simplicity  and  consistency  with  observations,  we  as¬ 
sume  the  rubble  accumulation  to  have  uniform  thick¬ 
ness  and  porosity,  and  consider  a  range  of  unit  ice 
volumes  t^l -ea)  from  0.75  m  to  1.5  m,  for  an  origi¬ 
nal  sheet  thickness  of  0.5  m. 

Ice  rubble  was  not  observed  near  site  D  during  or 
after  the  ice  motion,  except  for  minor  accumulations 
along  the  banks.  Prior  to  700  s  the  moving  ice  at  site 
U  was  also  large  plates  ora  continuous  sheet  with 
an  undisturbed  unit  volume  of  05  m.  The  quantity 
of  ice  accumulating  in  the  reach  between  the  initial 
motion  and  460  s  was  obtained  from  eq  15.  We  as¬ 
sume  that  half  of  this  total  resides  in  the  primary  ac¬ 
cumulation,  and  calculate  its  initial  length  L^,  from 
eq  17  as  a  function  of  unit  ice  volume.  All  of  the  ice 
entering  the  reach  after  460  s  is  assumed  to  be  part 
of  this  accumulation.  Between  700  and  860  s  the  ice 
passing  site  U  was  more  broken  and  concentrated, 
and  we  estimate  the  unit  volume  at  0.66  m.  At  860  s 
the  convergence  front  arrived  at  site  U,  and  the  unit 
volume  was  that  of  the  accumulation.  With  these  data 
the  final  accumulation  length  Lr  can  be  found  from 
eq  17  for  each  assumed  unit  ice  volume.  The  param¬ 
eters  Rb  and  C^/Vx  are  obtained  for  each  unit  vol¬ 
ume  from  eq  3  and  4,  and  all  these  results  are  sum¬ 
marized  in  Table  2.  Th<_  linensionless  breaking  front 
speed,  and  the  length  and  volume  of  the  rubble  vary 
significantly  with  the  assumed  unit  ice  volume. 

Three  hours  after  initial  ice  motion  a  second  sup¬ 
port  breakup  was  recorded  by  the  video  camera  at 
site  D.  This  videotape  indicated  that  most  of  the  ice 
rubble  in  the  reach  was  in  the  primary  accumulation, 
and  was  the  basis  for  that  assumption.  Ice  velocity 
data  were  obtained  and  polynomials  were  fitted  to 
these  data.  Integration  of  the  polynomials  provided 
length  scales  that  were  used  to  obtain  an  estimate  of 
the  length  of  the  earlier  ice  accumulation  between 
the  camera  sites.  The  river  at  the  dr"  vnstream  site  is 
narrower  than  in  the  reach  between  the  sites.  Cor¬ 


recting  for  the  difference  in  width  we  obtain  an  ini¬ 
tial  accumulation  front  location  1110m  upstream  and 
an  initial  accumulation  length  of  about  490  m.  Table 
2  indicates  close  agreement  of  this  length  with  that 
for  a  unit  ice  volume  of  1.10  m.  During  this  second 
ice  release  about  1  km  of  ice  sheet  and  8  km  of  rubble 
moved  by  the  downstream  site  at  an  average  speed 
ofalmost3m/s.  Withthel.l  accumulation  unit  vol¬ 
ume  used  in  eq  8  we  obtain  an  average  ice  discharge 
of  about  450  mVs. 


KINEMATIC  ANALYSIS  OF 
CONNECTICUT  RIVER  DATA 

Event  sequence 

We  will  now  develop  an  event  sequence  from  the 
data  with  the  ice  continuity  equations  in  order  to 
learn  as  much  as  possible  about  the  dynamics  of  the 
1992  Connecticut  River  breakup.  Throughout  the  ini¬ 
tial  420  s  of  motion,  the  ice  moved  as  large  plates  or 
sheets  at  both  sites  (Fig.  8, 9).  The  arrest  of  the  mo¬ 
tion  at  site  D  at  420  s  initiated  a  stoppage  front  that 
traveled  upstream.  This  front  arrived  555  s  later  at 
site  U,  corresponding  to  an  average  speed  of  only 
2.9  m/s.  With  this  front  speed  and  an  average  ice 
velocity  of  about  1  m/s  (Fig.  10),  the  downstream 
ice  must  initially  have  a  much  greater  unit  volume 
than  that  upstream  (Fig.  5).  However,  all  ice  visible 
from  site  D  was  a  uniform  sheet,  indicating  that  a 
single  stoppage  front  is  not  consistent  with  the  data 
and  observations.  The  abrupt  increase  in  ice  velocity 
at  site  U  after  460  s  corresponds  to  a  rapid  decrease 
in  the  local  resistance  to  motion.  The  initiation  of  a 
breaking  front  downstream  would  generate  this  re¬ 
sponse.  Therefore,  we  will  assume  that  the  stoppage 
front  traveled  rapidly  upstream  to  a  location  near  site 
U  in  abc  ut  460  -  420  =  40  s.  The  average  front  speed 
of  35  m/s  requires  negligibly  thicker  ice  downstream 
and  minimal  convergence  would  have  occurred.  At 
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460  s  the  stoppage  front  stalled  and  then  began  mov¬ 
ing  downstream  as  a  second  breaking  front.  The  ice 
accumulation  upstream  of  the  breaking  front  was 
thicker  and  had  a  greater  unit  volume  than  the  sta¬ 
tionary  plates  downstream.  The  convergence  front 
arrived  at  site  U  from  downstream  at  860  s.  Prior  to 
that  time  die  ice  velocity  decreased  almost  linearly 
without  noticeable  convergence.  If  site  U  was  typi¬ 
cal  of  nearby  locations,  the  steady  unit  ice  volume 
indicates  rapid,  rigid-body  response  to  changes  in  the 
motion  downstream.  We  hypothesize  that  after  460 
s  the  front  and  ice  motions  in  the  reach  downstream 
are  reflected  in  the  observations  and  velocity 
measurements  at  site  U,  and  that  the  regions  of  ice 
convergence  are  limited  to  the  fronts.  For  these  con¬ 
ditions  the  kinematic  model  indicates  that  the  break¬ 
ing  front  moved  continuously,  stalled  briefly,  moved 
again,  and  finally  arrested.  The  ice  in  the  ac¬ 
cumulation  mimicked  this  behavior.  While  the  break¬ 
ing  front  was  moving  downstream  the  convergence 
front  moved  upstream,  and  the  accumulation  length¬ 
ened  rapidly.  The  convergence  front  gradually  be¬ 
came  a  stoppage  front  as  the  breaking  front  speed 
and  accumulation  ice  velocity  decreased. 

Ice  and  front  motion 

We  will  now  calculate  the  ice  and  front  motion 
downstream  of  site  U  from  460  s  until  the  motion 
arrests  at  975  s.  For  simplicity  the  piecewise  linear 
approximation  of  the  ice  velocity  shown  in  Figure 
10  will  be  taken  as  Viu.  The  quantity  of  ice  obtained 
from  eq  15  and  the  linear  approximation  is  within 
1%  of  that  using  the  polynomial  for  a  wide  range  of 
accumulation  unit  volumes.  The  assumed  rigid  body 
motion  of  the  ice  between  site  U  and  that  just  up¬ 
stream  of  the  convergence  front  relates  Vic  to  Piu  by 
the  relative  river  widths 

Vic  =  ^3-  (19) 

B 

and  Vlc  =  of  the  CV  at  the  convergence  front.  The 
ice  accumulation  behind  the  breaking  front  is  also 
assumed  to  be  rigid  with  uniform  thickness,  poros¬ 
ity  and  speed,  and  the  ice  velocity  upstream  of  the 
breaking  front  Pib  =  Va  of  the  CV  at  the  convergence 
front.  The  rigid  body  assumptions  limit  the  ice  con¬ 
vergence  to  the  breaking  and  convergence  fronts. 

Measured  ice  velocities  on  either  side  of  the  con¬ 
vergence  front  together  with  known  unit  ice  volumes 
in  eq  5  yield  the  time  variable  front  speed.  However, 
here  we  have  only  a  single  velocity  measurement  and 
must  assume  an  average  convergence  front  speed  in 
order  to  calculate  the  ice  velocity  on  the  other  side  of 
the  front.  The  corresponding  ice  discharges  can  then 


be  obtained  with  eq  4,  9  and  11.  Parameter  ranges 
will  be  used  to  assess  the  sensitivity  of  the  calcula¬ 
tions  to  unit  ice  volume  of  the  accumulation  (0.83, 
1.1,  1.4  m)  and  initial  average  convergence  front 
speed  (-0.3,  -0.1,  0.1  m/s).  This  speed  is  taken  as  a 
constant  until  a  change  occurs  in  the  unit  ice  volume 
arriving  at  the  front  from  upstream.  We  assume  that 
this  change  should  not  affect  the  speed  of  the  break¬ 
ing  front.  The  ratio  of  the  ice  velocities  on  either  side 
of  the  convergence  front  is  set  equal  to  that  from  the 
previous  time  step,  and  Cc  is  recomputed  using  eq 
5.  This  assumption  can  be  visualized  in  Figure  3  as  a 
change  in  the  ratio  of  unit  ice  volumes  along  a  line 
of  constant  Vlb/  Vic,  providing  a  new  Cc/  Vic.  After¬ 
ward,  Cc  is  again  held  constant  and  ice  velocity  is 
computed  until  the  next  change  in  unit  volume  of 
incoming  ice.  With  the  time  of  convergence  front  ar¬ 
rival  at  site  U,  Cc  and  Lro  known,  the  initial  breaking 
and  convergence  front  positions  can  be  determined. 
Then,  the  front  positions  and  all  ice  particle  positions 
between  sites  U  and  D  can  be  determined  through 
time. 

Ice  velocity  and  ice  discharge  upstream  of  the 
breaking  and  convergence  fronts,  and  celerity  of  the 
breaking  front  are  presented  in  Figure  12  for  ranges 
of  accumulation  unit  volume  and  initial  convergence 
front  speed.  The  zero  time  in  Figure  12  corresponds 
to  an  event  time  of  460  s.  The  ice  accumulation  ve¬ 
locity  V,b  increases  as  its  unit  volume  decreases  and 
as  the  initial  front  speed  increases.  Vlb  is  obtained 
directly  from  the  measurements  after  400  s  (860  s 
event  time),  and  does  not  vary  with  these  parameters. 
The  ice  velocity  upstream  of  the  convergence  front 
V'to  obtained  from  eq  19  and  the  measured  data,  was 
always  greater  than  Vjb  prior  to  400  s  and  unknown 
afterward.  The  ratio  Vib/V ic  <  1  indicates  ice  conver¬ 
gence  at  the  upstream  end  of  the  accumulation  over 
the  assumed  parameter  ranges.  The  trends  in  break¬ 
ing  front  celerity  follow  those  of  the  ice  velocity.  The 
exception  is  after  400  s  when  Cb  increases  with  de¬ 
creasing  accumulation  unit  volume,  while  Pib  does 
not  have  this  dependence.  The  ice  discharge  of  the 
accumulation  upstream  of  the  breaking  front  Qt,  in¬ 
creases  slightly  as  the  accumulation  unit  volume  de¬ 
creases,  and  increases  more  significantly  as  the  ini¬ 
tial  convergence  front  speed  increases.  After  400  s 
Qt,  varies  only  with  the  unit  volume  and  displays  the 
opposite  trend.  The  ice  discharge  upstream  of  the 
convergence  front  Qc  does  not  vary  over  the  param¬ 
eter  ranges.  Differences  in  Qc  result  from  the  vari¬ 
able  arrival  times  of  thicker  ice  from  upstream  at  the 
convergence  front.  Qc  is  always  greater  than  Qt,  ex¬ 
cept  when  Cc  is  positive. 

A  comparison  of  front  speeds  and  positions  is 
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Figure  12.  Ice  velocity  and  ice  discharge  upstream  of  the  breaking  and  convergence  fronts,  and  the  celerity  of  the  break¬ 
ing  front,  for  a  range  of accumulation  unit  volumes  and  initial  convergence  front  speeds. 


given  in  Table  3  for  the  assumed  parameter  ranges. 
The  initial  convergence  and  breaking  front  positions 
Xco,  Xbo  were  always  near  site  U  at  x  =  0  with  a  total 
range  of  150  m.  The  intersection  time  of  the  thicker 
ice  from  upstream  with  the  convergence  front  in¬ 
creased  with  decreasing  convergence  front  speed  and 
accumulation  unit  volume,  having  a  total  range  of 
50  s.  The  change  in  convergence  front  speed  follow¬ 


ing  the  intersection  A Cc  decreased  with  increasing 
unit  volume,  but  was  insensitive  to  initial  conver¬ 
gence  front  speed.  The  total  distance  traveled  by  the 
breaking  front  AX*>  and  the  average  breaking  front 
speed  Cb  increased  with  convergence  front  speed  and 
with  decreasing  unit  volume  of  the  accumulation. 
The  positions  of  the  breaking  and  convergence  fronts 
and  several  ice  particles  are  traced  through  time  on 
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Table  3:  Comparison  of  front  speeds  and  positions  for  ranges  of  accumulation  unit  volume 
and  initial  convergence  froc*  speed. 
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Figure  13.  Breaking  and  convergence  front  and  ice  particle  motion 
on  the  x-t  plane  that  are  most  consistent  with  the  Connecticut  River 
data. 


the  x-t  plane  in  Figure  13.  The  ice  motion  depicted 
is  that  most  consistent  with  the  data,  Cco  =  -0.1  m/s 
and  fir  (1  -  e^)  =  1.1  m.  The  ice  particle  that  departs 
from  site  U  at  240  s  (event  time  700  s)  represents  the 
thicker  ice  from  upstream  that  produces  a  kink  in 
the  convergence  front  trace.  As  observed,  the  com¬ 
puted  position  of  the  convergence  front  passes  site 
U  moving  upstream  at  400  s.  Figures  4  and  13  both 
depict  developing  ice  accumulations.  However, 
changing  ice  and  front  speeds  produce  paths  with 
variable  slopes  that  replace  the  parallel  lines  of  steady 
motion. 


CONCLUSIONS 

We  developed  a  kinematic  model  of  ice 
motion  during  dynamic  river  breakup  by 
applying  continuity  at  several  fronts  in  the 
ice.  Breaking,  convergence,  stoppage  and 
release  fronts  were  described,  and  relation¬ 
ships  were  developed  between  their  respec¬ 
tive  speeds  and  the  changes  in  unit  ice  vol¬ 
ume  and  ice  discharge  across  each  front. 
Rapid  breaking,  stoppage  and  release  front 
propagation  indicate  minimal  modification 
of  the  ice  field,  while  lower  speeds  corre¬ 
spond  to  larger  changes.  The  difference  in 
ice  speed  and  unit  volume  across  the 
convergence  front  determine  the  front 
speed  and  direction.  A  convergence  front 
becomes  a  stoppage  front  when  the  ice  ac¬ 
cumulation  speed  is  zero.  Extended  break¬ 
ing  front  motion  at  speeds  comparable  to 
the  ice  velocity  causes  the  development  of 
an  accumulation  while  the  ice  is  in  motion. 

The  kinematic  model  was  applied  to  a 
dynamic  breakup  of  the  Connecticut  River 
using  measured  ice  velocity  data  at  sites 
both  upstream  and  downstream  of  a  study 
reach.  Important  assumptions  were  that 
changes  in  ice  and  motion  conditions  oc¬ 
curred  only  across  a  front,  and  that  a  developing 
accumulation  would  have  uniform  thickness  and 
porosity.  The  kinematic  analysis  supports  the  rela¬ 
tionship  described  by  Ferrick  et  al.  (1992)  between 
ice  convergence  upstream  of  a  breaking  front  and 
front  speed.  The  extremely  fast  measured  breaking 
front  speed  that  initiated  the  ice  motion  in  the  reach 
supports  the  force  balance  per  unit  length  and  the 
instantaneous  support  breakup  of  a  river  reach  in  the 
model  of  Ferrick  and  Mulherin  (1989).  The  analysis 
together  with  the  data  provided  a  hypothetical  event 
sequence  of  ice  processes  in  the  reach.  The  accumu- 
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lation  length  and  volume  were  obtained  as  functions 
of  the  unit  accumulation  volume,  and  the  mean  ac¬ 
cumulation  unit  volume  was  estimated.  Ice  and  front 
motions  and  ice  discharge  were  determined  through¬ 
out  the  reach  for  ranges  of  accumulation  unit  vol¬ 
ume  and  initial  convergence  front  speed.  The  mo¬ 
tion  simulation  most  consistent  with  all  the  data  and 
observations  was  presented  in  the  x-t  plane.  Accord¬ 
ing  to  this  analysis,  the  ice  accumulation  formed 
while  in  motion,  and  then  arrested,  creating  a  jam. 

This  kinematic  model  of  ice  motion  during  river 
breakup  represents  a  new  approach  to  the  develop¬ 
ment  of  a  dynamic  model  that  includes  jam  forma¬ 
tion  and  release.  Static  models  of  ice  jam  formation 
neglect  both  dynamic  ice  motion  and  unsteady  river 
flow.  In  contrast  to  the  static  models  that  are  essen¬ 
tially  independent  of  breakup,  the  kinematic  model 
details  the  breakup  processes.  Static  models  are  only 
applicable  when  the  speed  and  duration  of  breaking 
front  motion  are  small  prior  to  arrest,  so  that  the  ef¬ 
fect  of  this  motion  on  the  force  balance  is  negligible. 
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